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The palmprint recognition system has gained significant attention in security 
and law enforcement due to its unique features, such as principle lines, 
ridges, and wrinkles. However, many existing methods for extracting these 
features have limited accuracy, especially when the image illumination 
varies or the size of the processed pixels increases. Previous studies have 
shown that the local binary patterns (LBP) algorithm is effective for 
palmprint recognition due to the rich texture characteristics of a palmprint. 
In this paper, we propose a new technique for a robust contact-based 
palmprint identification system using vertical-LBP and KAZE feature 
detection. Our technique aims to improve recognition accuracy by using 
KAZE, which is a nonlinear diffusion approach that extracts nonlinear 
features from the evolution of the illuminance of an image. We also utilize 
principal component analysis (PCA) to reduce the dimensionality of the 
generated descriptor vector elements. The proposed method was tested on 
the PolyU database and achieved recognition accuracy of 99.7%. 
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1. INTRODUCTION 


Palmprint recognition is considered the most efficient technique of biometric recognition, with the 
advantages of low cost, low-resolution, non-intrusion, reliability, and high performance in the field of 
biometric recognition [1], [2]. There are two types of palmprint images: high-resolution and low-resolution. 
Each type is appropriate for specific applications; for example, forensic applications usually use 
high-resolution images, while low-resolution images are used for access controls [3]. 

Since Shu and Zhang [4] first proposed it, palmprint has piqued the interest of numerous researchers 
numerous methodologies have been introduced and expanded for palmprint recognition [5] which are 
primarily based on: texture, principal lines, and subspace. The proposed method primarily focuses on 
designing texture-based local binary pattern (LBP) due to its reliability. LBP which has been proposed by 
Ojala et al. [6] gathers momentum in the research area of palmprint recognition and several LBP variants 
have been proposed. These variants include local ternary pattern (LTP), dominant local binary pattern 
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(DLBP), center-symmetric local binary pattern (CS-LBP), local derivative pattern (LDP), and completed 
local binary pattern (CLBP) [5]. 

Sehgal [7] proposed a palmprint recognition system using LBP and support vector machines (SVM), 
while Tamrakar and Khanna [8] presented an approach to palmprint identification based on local uniform 
LBP (ULBP) distribution. ULBP applies to the pattern with minimal discontinuities of uniform appearance. 
Tarhouni et al. [9] proposed a feature extraction technique using a pyramid-oriented gradient histogram with 
LBP, where the characteristics for the classification of palmprints are concatenated, while multi-resolution 
capabilities could represent more image features in different settings. Guo et al. [10] proposed a collaborative 
representation model with hierarchical multi-scale LBP (HMS-LBP) for palmprint recognition. HMS-LBP 
can obtain useful information from non-uniform patterns as well as from the grayscale, rotation, and 
illumination effects. Dubey and Kanumuri [11] an orientation diagram was extracted using anisotropic filters 
(AF) from an optimal local binary directional pattern (OLdirBP). The characteristics were finally achieved 
using local directional patterns accompanied by LBP. Dubey et al. [12] proposed the ability of the binary 
wavelet transform (BWT) to discriminate against palmprints. Further, micro and macro-pattern histograms 
are extracted using local binary pattern (LBP) from different transformed bit planes and concatenated to form 
the feature vector. Several recent studies [13], [14] have suggested combining hierarchical multi-scale 
complete local binary pattern (HMS-CLBP) and weighted sparse representation-based classification 
(WSRC). However, basic LBP it is not rotation invariant, meaning that the LBP features will change if the 
image is rotated. This can be addressed by using extended LBP or by using vertical local binary patterns 
(VR-LBP), in combination with other descriptors that are rotation invariant. 

Attallah [15] proposed a novel approach for the palmprint fusion of multiple features. These features 
are provided from multi-spectral images at 940 nm, which allow the extraction of the information under the 
skin of the palm. In this context, CLBP, which added an extra bit for each bit encoded by LBP corresponding 
to a neighbor of the local neighborhood, was used to construct a robust feature descriptor that exploits both 
the sign and the inclination information of the differences between the center and the neighbor grey values. 

There are many methods for local feature detection, and the most common are speeded-up robust 
features (SURF) [16], scale-invariant feature transform (SIFT) [17], and oriented fast and rotated brief (ORB) 
[18]. They detect and define features at different scale levels by constructing or approximating the image's 
Gaussian scale space. Moreover, the Gaussian blur violates the natural limits of images and smooths all 
information with noise to the same point, decreasing the precision and distinctiveness of the localization [19]. 

Alcantarilla et al. [20] introduced the KAZE feature, which is a wind-meaning to extract features with 
higher discrimination power; the wind is characterized in nature as a large-scale flow of air, and thus this flow is 
usually regulated by nonlinear processes. The KAZE feature detects and describes 2D features in nonlinear 
scale space using nonlinear diffusion filtering. The KAZE adjusted blurring locally to the images, reducing 
noise while maintaining object borders, achieving superior precision and distinct positioning. Figure 1 illustrates 
the four steps of the classic method for palmprint recognition: getting an image, prepossessing it, feature 
extraction, and matching. 


Figure 1. Palmprint recognition system 


The main contribution of this paper is proposing a new technique that combines VR-LBP and 
KAZE features, which further improves the accuracy of recognition by generating higher knowledge of 
discrimination during the extraction of the features by using PolyU database contact palmprint images. 
Among the many emerging techniques for palmprint recognition, the high recognition levels produced by 
various palmprint image texture feature extraction methods are the most common. The following is the 
structure of this paper: section 2 describes the method, section 3 highlights the results and discussion of the 
experiment, and section 4 is the conclusions. 


2. METHOD 

The palm has texture features, and to extract these features, the algorithm must be modified to 
extract more features from the palm. In the proposed method, the way to compare the neighbors with the 
center pixel of VR-LBP is changed to a vertical one rather than the default one; however, the grayscale 
palmprint image is typically a texture image that primarily consists of ridges and principal lines. By using 
VR-LBP, it will ensure recognition accuracy by having many more features than the basic LBP, as shown in 
Figure 2. 
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Figure 2. The basic operator of LBP 


Here is the difference between VR-LBP and the original one: by comparing the central pixel with its 
neighbors in a vertical way, the value of the central pixel will be different, as shown in Figure 3, where 0-7 are 
the sequences of neighbors in a vertical way using (1). 


LBPp.p= D5-0 S(9p — Gc)2” (1) 
s(x) = ka 2) 


Where g. is gray value of a central pixel, gp is the value of neighbors, P is the total number of involved neighbors, 
R is the radius of the neighbors, and gp is (Rcos(27p/P), Rsin(27p/P), Figure 4 shows the basic LBP vs VR-LBP, 
where Figure 4(a) is the original palm image, Figure 4(b) applying LBP and Figure 4(c) applying VR-LBP. 
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VR-LBP binary code 


Figure 4. The different between LBP vs VR-LBP; (a) original palm image, (b) basic LBP, and (c) VR-LBP 


To promote recognition accuracy, our method is applied, and it can be summarised in four steps, which 
are: i) the final function code is derived from VR-LBP, while the number of neighbors for VR-LBP that has 
been used is 3x3, ii) extract the features by using KAZE feature detection in different parameters, iii) apply 
PCA for dimension reduction, which is proposed in [21], and iv) finally, the produced features are matched 
using euclidean distance. Figure 5 depicts the proposed method's flow chart, and Figure 6 depicts an example of 
using VR-LBP. 


Bulletin of Electr Eng & Inf, Vol. 13, No. 2, April 2024: 1060-1068 


Bulletin of Electr Eng & Inf ISSN: 2302-9285 O 1063 


PCA 
Dimension reduction 


MATCHING 
Euclidean distance 


Figure 5. Proposed method flow 
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Figure 6. Samples of VR-LBP 


Following the implementation of VR-LBP, the following steps are taken to detect the KAZE feature: 
first, build an additive operator splitting (AOS) input image non-linear scale space and variable conductivity 
diffusion system. Second, the calculation of the scale-normalized Hessian matrix in the nonlinear scale space. 
Third, detect the 2D feature points in the special Hessian matrix by measuring the local limit of 3x3 
neighborhoods; last but not least, calculate the feature points and obtain an invariant scale and rotation 
descriptor dependent on the image derivatives of the first order [22], [23] nonlinear diffusion as (3): 


< = div(c(xy,t).L) (3) 


where V and div are respectively the divergence and gradient operators and c is the conductivity function of 
the diffusion equation, time (t) is the scale parameter, and larger values lead to simpler image representations. 


c(xy,t) = g(WVLa(xy,t)|) (4) 
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The VLo is the gradient of a gaussian smoothed version of the original image L. Perona and Malik described 
two different formulations for the conductivity function g: 


_ [WLal? [VLal? 


g1 = exp ( Meat) g2 = 1/1 + aL (5) 


where k is the contrast factor regulating the diffusion level. The g1 feature encourages edges of high contrast, 
while g2 encourages large and smaller regions [20]. Additive operator splitting (AOS) formulation: 


pit1_yi 


= Xk, A Cn (6) 


T 


where A; is a matrix encoding each dimension's image conductivity. The solution L! can be obtained as (7): 
Lt = (1 -T YE AL) (7) 
For detecting points of interest: 
LyHessian = a? (LyxLyy SS L) (8) 


where (Lxx, Lyy) are the horizontal and vertical derivatives of the second order, and where Lyy is the cross 
derivative of the second order. Due to the set of filtered images from Li nonlinear space. The keypoints 
location is determined using the proposed in [23] sub-pixel accuracy. 

For the dimension reduction using linear projection based on the Karhunen—Loeve approach, the 
PCA technique reduces the dimensionality of palmprint data. This approach maximizes the distribution of all 
the predicted image samples of palmprint. Dimensional feature reduction increases the power of 
discrimination and reduces the time and memory cost of computing. In the PCA algorithm, the linear 
projection approach is used to project high-dimensional data into a low-dimensional subspace [24]. 

Consider a collection of sample palmprint images {y1, y2, yn} collected image space to illustrate the 
computation of main components, and assume that each palmprint image belongs to a different class of 
classes {X1, X2, Xc}. Consider the linear transformation that transforms the original high-dimensional palmprint 
image space with n-dimensions into a new m-dimensional low-dimensional feature space where m<<n 
obtained from (9) which is the linear transformation equation, may represent the generated feature vectors y“=R™. 


Yr = Wx (9) 
where k=1,2, N, (10) represents the entire Scatter matrix Sr. 
Sp = DVRe1 Oe — De — W)" (10) 


where N is the total sample of palmprint images and u E R™which is the mean vector of all palmprint 
images, W'SrW is the output of yx by linear transition with WT. T optimizing the scatter matrix determinant 
obtained from the projected image of the palmprint, and projection Wop is chosen to maximize the 
determinant of the expected sample total scatter matrix, i.e.: 


Wopt = argmax|WTS rW] = [w1Ww2....Wm] (11) 


where [w1 W2...Wm] represents the eigenvectors with n-dimensions from which they are chosen based on their 
own corresponding largest values. 


Z(m) =T(m) —X (12) 
where Z is the 2D zero-mean training data matrix, (X) is the mean row vector of the training data matrix, (m) 


is the average palmprint image vector of a given subject's training palmprint, and T is the amount of the 
variable. For eigenvalues and eigenvectors from the covariance matrix using as (13): 


S = EE" (13) 


where S is the 2D covariance matrix, x is the zero-mean 2D training data matrix, and xT is the zero-mean 2D 
training data matrix transposing matrix. Formula for measuring the eigenvectors and eigenvalues is as (14): 


S.e = Xe (14) 
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where S is the 2D matrix, e is the eigenvector matrix, and A is diagonal matrix of the eigenvalue. 


eS Z* projected (15) 
Where projected is the eigenvector after projection and F is the feature matrix of training data, Z is the zero 
mean data matrix. At the last step for matching, the euclidean distance classifier, to calculates the distance 
between features in the template and recognition according to the smallest distances. In a cartesian coordinate 
system, assume that g=(q1, q2, qn) and p=(p1, p2, Pn) are two points in the euclidean n-space. Then, by using 
the Pythagorean (16) [25], the example of the euclidean distance between two points is shown in Figure 7. 


ED(p,q) = VÈ; (pi — qi)? (16) 


Testing Samples Training Template 


Figure 7. Euclidean distance 


3. RESULTS AND DISCUSSION 

In this section, we have developed a new technique for a robust palmprint identification system 
using VR-LBP and KAZE feature detection. The PolyU database has been used to validate our proposed 
method’s efficiency and compare it with other palmprint recognition algorithms. The database includes 7,752 
palmprint images with a size of 384*286 from 386 palms. Palmprint images have been collated in two 
different sessions, containing more than nine images. Two sessions have an average duration of two months, 
introducing improvements in light conditions, emphasis on the imaging system, and skin texture deformation 
due to various hand positions between the two image capture occasions. 

Assumedly, accuracy is one of the most important considerations to be considered. Thus, in this 
segment, most analyses concentrate on accuracy rather than other factors, even though, during the 
comparison, our methods were compared with those of other works. In these experiments, an evaluation 
protocol has been used, like nearest feature detection according to our detection method (KAZE), and 
compared with SURF based on the number of feature points as shown in Table 1. From the Table 1, KAZE 
feature detection is better than SURF in terms of detection points at the default parameters, which can help 
ensure recognition accuracy according to the increasing number of features. 

In this work, the sharp edge was chosen according to the number of points that can be seen in Table 2. 
Nevertheless, in the VR-LBP, the number of points is increasing, which means it's improved and will improve 
the accuracy during the matching process. There are several parameters to be defined in the proposed algorithm 
that will affect the recognition rate, such as the diffusion value, the sub-image dimension, the value of A, and the 
PCA dimensions. In this part of the experiment, all classes (386 palms) are used; for each palm, 8 samples are 
used as training samples, and | sample for each class is used as a testing sample. 


Table 1. Comparison of KAZE and SURF Table 2. Comparison KAZE with changing the value of diffusion 


"E KAZE KAZE Original image VR-LBP Basic LBP 
SURF (points) (points) parameters (points) (points) (points) 
Original image 0 610 Edge 198 631 628 
VR-LBP 276 634 Region 610 634 628 
Basic LBP 257 628 Sharp-edge 203 657 632 


In Table 3, could notice the accuracy increasing while the number of classes and the number of 
images increased according to the PCA algorithm and with the improvement of VR-LBP+KAZE. From 
Table 4, using basic LBP and PCA, the accuracy is less without KAZE compared to our proposed method, 
and from Table 5, it appears that the method of Zhang [1] and El-Tarhouni [9], from Table 4, using basic 
LBP and PCA, the accuracy is less without KAZE compared to our proposed method, and from Table 5, it 
appears that the method of Zhang [1] and El-Tarhouni [9], is less accurate since the number of classes that 
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have been used is not that much the same as our proposed method; nevertheless, increasing the number of 
neighbors will lead to an increase in time and computation complexity. In the Tables 3 and 4 can see the 
accuracy increasing as the number of classes and images increases. This is due to the PCA algorithm and the 
improvement of VR-LBP+KAZE. The fact that the accuracy is lower than that of VR-LBP demonstrates that 
the proposed method is effective. Figure 8 depicts the accuracy of the basic LBP and VR-LBP when KAZE 
and PCA are used. The experiment was carried out with MATLAB 2018b installed on a computer with a 
Windows 10 platform, a CPU core i178565U, 8 GB of RAM, and 533 MB/s read and 442 MB/s write speed 
solid-state drive (SSD) memory. 


Table 3. VR- LBP with KAZE sharp-edge diffusion Table 4. Basic LBP+KAZE + PCA 
+ PCA using euclidean distance 
Class number Number of train image Accuracy (%) Class number Number of train image Accuracy (%) 
386 2 76.16 386 2 66.7 
3 84.19 3 72.9 
4 90.67 4 83.5 
5 94.81 5 87.9 
6 97.68 6 91.34 
7 98.44 7 91.6 
8 99.7 8 93.6 
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Figure 8. Accuracy graph between two methods 


Table 5. Comparison of previous work 
Method Number of class _ Accuracy (%) 


[1] 129 98.06 
[9] 250 99.27 
Our proposed 386 99.7 


4. CONCLUSION 

To avoid computation complexity, the discrimination capability of the image is used in this work for a 
robust palmprint identification system using VR-LBP with several neighbors (3*3). The system performs faster 
compared with the different number of neighbors, VR-LBP which made the palm features more apparent than 
basic LBP which leads to an increase in the recognition accuracy, nevertheless, KAZE feature detection used 
nonlinear diffusion approaches for detecting the palm features which confirmed the result when compared with 
other methods, last but not least use PCA for reducing the dimensionality of the generated descriptor vector 
elements. The proposed method was tested on the PolyU database and successfully improved the recognition 
accuracy, resulting in a high percentage of 99.74%. 
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